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Оцщнка похибки заокруглення наближення функщй 
двох зм1нних одном!рними операторами 


В статье предложен подход к оценке ошибки округления, возникающей при приближении функции двух 
переменных одномерными операторами смешаннои аппроксимации при решении системы линеиных 
алгебраических уравнений с помощью обратной матрицы. 
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В статт! запропоновано шджд до ощнки похибки заокруглення, що виникае при наближенн! функцИ двох 
зманних одновимрними операторами, шляхом розв’язання системи лиййних алгебрайчних равнянь, що 
зводиться до обчислення обернено! матрицу. 

Ключов! слова: одновимрний оператор, змппана апроксимащя, похибка заокруглення, 
похибка методу, похибка даних, обернена матриця 


Введение 


С каждым годом сложность проведения вычислительных расчетов на ЭВМ 
увеличивается, возникает необходимость решения научных и прикладных задач, 
связанных с обработкой больших массивов входных данных. Все более актуальной 
становится проблема оценки ошибки округления при приближении функции двух и 
более переменных одномерными операторами, возникающей при расчетах на ЭВМ, 
и определение ее влияния на общий результат с учетом ошибок метода вычисления и 
входящих данных. 
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Сегодня широко используются интерполяционные и аппроксимационные методы 
приближения функций, среди которых можно отметить приближение функций двух 
переменных, заданных следами на системе взаимно перпендикулярных прямых или 
набором дискретных данных через операторы, действующие на функцию по одной 
переменной [1]. В дальнейшем будем называть такие операторы одномерными. 

Цель работы — разработать метод оценки ошибки округления приближения 
функции двух переменных заданных следами на системе взаимно перпендикулярных 
прямых, через одномерные операторы, рассмотренные в работе [1], учитывая вид 
операторов и предположение, что соответствующие интегралы вычислены точно. 
Считаем, что ошибка округления возникает в ходе решения систем В\и = ЕТ (у) и 


уВь = Е› (х) путем нахождения решения с помощью обратных матриц. 


Постановка задачи 
В работе [2] авторы доказали теорему об ошибке приближения, включая ошибку 
метода приближения и ошибку входящих данных, функции двух переменных 


У (х, у) е С? [0/1]? оператором сплайн-интерлинации вида 


И м м мм _ 
ОГ(х, у) = УФ + Уи (0) = Ут), 
&=0 [=0 


к=0/=0 
когда следы /(х,,у), Г(х,У,) заданы функциями ф; (у), и; (х) с ошибками и зна- 
чения /(х,,у/) также заданы числами ДЕ с ошибками. Рассмотрим подходы к 


оценке ошибки округления при применении разработанного в [2] метода приближе- 
ния функций двух переменных, заданных следами на системе взаимно перпендику- 
лярных прямых, когда оператор (1) задан в виде 


2 Л(х, у) = АУ (х, у) + Л, у) - АА Л(х, У), (2) 
где одномерные операторы 4А|, 45 задают наилучшее среднеквадратическое прибли- 
жение функции /(х, у) по переменным х и у, и определяются формулами 


А (ху) = И (х)ф' (р) = ПОВЕИЕГ (), А Г(х, у) = (хо (У) = Е, (В ©, 
А! А Х(х, у) = В (ВЕ ЕВ Ва (У), вн (х) = [о (0)... МОЛЬ 


1 1 
5 (7) = [15.0 СУЛЬ В, = [м о дах, Вь = [№ (№5 ау, 
0 0 
1 1 
Е= [м (ГУ) (оахау, ВГ (= [Г 9. сь дах, В (х) = [ (к, УР (у)ау. 
[в 0 0 


Теоретические подходы к оценке ошибки округления 


Отметим основополагающие труды, посвященные оценке ошибки округления [3-8]. 
В работе [9] предлагается метод оценки погрешностей округления, отличный от рас- 
смотренных в [3-8]. 

Ошибка округления — это ошибка, возникающая при выполнении арифметиче- 
ских операций на ЭВМ с округлением результатов до фиксированного количества 
разрядов. Различают два режима работы ЭВМ - с фиксированной запятой и пла- 
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вающей запятой. Для расчетов с фиксированной запятой каждое число находится в 
интервале —1<х<1, в который исходные числа приводятся путем масштабирования. 


При расчетах с плавающей запятой каждое число представляется в виде х = 22а ‚гдер — 
целое положительное или отрицательное число, называемое порядком, а (мантисса) — 


число, удовлетворяющее одному из неравенств: —-1<а<- и ИЛИ р <а<1 [3]. 


Теоретические сведения по оценке погрешности округления описаны в [5], [6]. 

В работе [3] комплексно анализируются важнейшие вычислительные аспекты 
определения ошибки математической модели и построения ее оптимальных реализаций. 
Комплексный подход основан на анализе трех основных характеристик вычислитель- 
ных методов — точности, времени реализации, требуемой памяти ЭВМ. По данным 
характеристикам выполняется сравнительный анализ и оптимизация соответствую- 
щих численных методов. Отмечено, что в практике численного вычисления задач на 
ЭВМ применяются следующие характеристики задач, алгоритмов и ЭВМ: Е(1,Х,У) — 
полная ошибка решения Е задачи Р на ЭВМ с помощью алгоритма А, Т(Г,Х,У) — 
время, необходимое для получения решения задачи; М (1, Х,У) — необходимая память 
ЭВМ; /е/ - коэффициент технико-экономической эффективности. 

В свою очередь полная абсолютная ошибка решения задачи Р(Г) на ЭВМ 


С(У) с помощью вычислительного алгоритма А(Х’) определяется так 


А(1,Х,У) = р(В, АУТ, = р(В,В,). 


А(Г,Х,У) = А +А) +Ду, Ао =р(Е А(Х)Т,), где А; = (К, Кол) — не 


=,4,й > 
устранимая ошибка решения задачи или ошибка за счет неточности входных данных; А _ 
Дз = р(К., А(Х,У)Г,) — ошибка 


округления реализации вычислительного алгоритма А(Х) наЭВМ С(Т). 


ошибка алгоритма А(Х) для определения К; ол, 


В дальнейших рассуждениях рассмотрим два способа нахождения обратной матри- 
цы, определяющие подход к оценке ошибки округления при приближении функции двух 
переменных оператором вида (2) с использованием одномерных операторов. 

В работе [7] отмечено, что значительная часть наиболее известных численных 
методов решения систем линейных алгебраических уравнений (далее — СЛАУ) 


Ах=Ь (3) 


основаны на разложении матрицы А на сомножители. 

В зависимости от того, как связаны сомножители с матрицей А, различают две 
схемы построения методов. В первой схеме предполагается, что явно известны сами 
сомножители, на которые разложена матрица А. Пусть 

А=В.С. (4) 

Решение системы (3) сводится к последовательному решению таких систем: 
Ву=Ь, Сх=у. 

Во второй схеме предполагается, что найденные матрицы [Ё,5,С, для которых 
выполняется соотношение 


[Аб =С. (5) 
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Тогда 
х=эи, (6) 
где и — решение системы 
Си =1 (7) 
с матрицею С из (5) и правой частью 
[=1Ь. (8) 


Решение системы (3) сводится теперь к вычислению вектора / согласно (8), 
решению системы (7) и определению искомого вектора х по формуле (6). В данной 
схеме матрицы Ё и 5 обычно бывают представлены в виде произведения элементарных 
матриц. 

Разложения (4), (5) можно использовать для вычисления обратной матрицы. 
Для (4) следует, что 


А = СВТ", (9) 


а из (5) имеем 


Ре (10) 
Поэтому если выполнено преобразование (4) или (5), то для получения матрицы 


1 
А’ остается только преобразовать одну или две матрицы простого вида и осуществить 
умножение матриц согласно (9) или (10). 
К задаче вычисления обратной матрицы можно подойти несколько иначе. Матрица 


—1 
А `` является единственным решением матричного уравнения АХ =Ё. 
Обозначим через х|,....х,„ вектор-столбцы матрицы А”'. Тогда х, является 


решением системы линейных алгебраических уравнений 
р (0 


где е; — координатный вектор с единицей на 1-м месте. Снова для решения системы 


(11) будут полезными разложения (4), (5). 

С практической точки зрения безразлично, вычислять обратную матрицу по форму- 
ле (9), (10) или с помощью решения систем (11). Автор работы [7] предпочитает второй 
способ, поскольку все вопросы, связанные с решением систем, уже исследованы. При 
реализации этого способа может потребоваться некоторое изменение вычислительной 
схемы методов, вызванное необходимостью одновременного решения систем (11) со 
многими правыми частями. 

Предположим, что выполнено преобразование (5), причем матрица С’ - правая 
треугольная, а матрица 5 представлена в виде произведения И\1...И„_›И„_ матриц 


отражения. Последовательное решение систем (7) с правыми частями е|.....е, вновь 


позволяет разместить всю информацию о решении на месте соответствующих столбцов 
матрицы С’. Преобразование (6) будем осуществлять последовательно, путем умно- 
жения сначала всех векторов на ЦИ„_|, затем на Ц„_› и, наконец, на И! . При этом 


следует учитывать как специальный вид преобразуемых векторов, так и специаль- 
м м —1 
ный вид самых преобразований. Элементы матрицы А’ снова могут быть получены 


—1 С 
на месте матрицы А ` после выполнения около Зи? арифметических операций. 
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Для численных методов решения системы линейных алгебраических уравнений, ос- 
нованных на разложении матрицы, для каждого столбца реально вычисленной матрицы 


— 


АГИИ для самой матрицы А”! можем записать [8] 


|" 
Е 1+1 
=. <2у.Л(тр о”, (12) 
| 
Е 
где у. — евклидовое число обусловленности матрицы А, у/ = | |“, 


—1 
Применение систем (11) для вычисления матрицы 2” позволяет, если будет 
необходимо, уточнить отдельные или все ее столбцы [8]. 
Кроме того, провести оценку ошибки округления построения обратной матрицы 


можно с помощью следующих соображений [7]. Пусть Т’- реально задана или реально 
вычисленная по некоторому вектору Б согласно предварительным условиям матри- 


ца вращения. Возьмем любой действительный вектор а. Если т и [а] <; Гораздо больше 


машинного нуля, где т -— количество базисных элементов, моделирующих число 
разрядов р -ичной системы исчисления вычислительной машины, то при умножении 


матрицы Г на вектор а выполняются соотношения 
Ла) =Та+ у || < Ура], Ла) = Раче) [в 2], 


где Л(.) - результат вычислений на ЭВМ выражения в скобках, {# — количество 


разрядов после запятой (в случае фиксированной запятой) или мантиссы (в случае 
плавающей запятой). 


Для матрицы Т, реально вычисленной, согласно тому же условию матрица 


ТТ является и и при этом выполняются следующие оценки 
1+1 


|77'-Е, > р - 15 МР = (22 +32)? =1+и, 


Пусть по ее Ь ее вычисляется матрица Т и вычисляется 
единственная ненулевая координата вектора ТЬ ‚ тогда 


::8 


И(ТЬ)=Т +2), РЕ. 


Приведенные оценки говорят о том, что т. вычисленная матрица вращения Г 
с высокой степенью точности не только близка к некоторой ортогональной матрице, 
но даже близка к ортогональной матрице, получаемой при точных вычислениях. При 
этом выявляются малыми и эквивалентные возмущения преобразованных векторов. 

Пусть и-мерный вектор 2 умножается на последовательность из № матриц 


— 


т, т Предположим, что Г, , > 5 „„/, - Реально заданные или реально 
вычисленные матрицы вращения, удовлетворяющие описанному выше условию. Тогда 


для любой последовательности пар индексов И Л.....йу/м имеют место соотношения 
—#+1 
ЗУМ". 


вращения Г, ГВ 


— 


ле, ое а у \ — 5), 
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Отметим, если точное решение есть результат реализации некоторого алгорит- 
ма над входными данными А, а приближенно вычисленное решение можно рас- 
сматривать как результата реализации того же точного алгоритма над входными 
данными А,, то отклонение А, от А называется эквивалентным возмущением [7]. 


Для эквивалентного возмущения М матрицы А при разложении А на множи- 
тели имеет место оценка 


АВЕО Т% И 


Тогда реально вычисленное решение х системы Ах = является точным ре- 
шением возмущенной системы (А-+&)х =Б +=. При этом 


В; Фр А, 


2: ур", 


где ф(и)+у(п) < 2Г(п), если только в пределах таких возмущений матрица остается 


невырожденной. 
При этом выполняется неравенство 


которое является следствием неравенства р: _ р и <21(тр "т |4] в 


ие <2у/(тр” (14) 


|8. 
Согласно (14) точность любого метода полностью определяется точностью раз- 
ложения матрицы на множители. 


Результаты вычислительного эксперимента 


Учитывая, что в операторе (2) используются одномерные операторы, дейст- 
вующие на одну переменную, проведем вычислительный эксперимент и определим 
ошибку округления при приближении функции одной переменной Х(х). 

На основе теоремы о виде остатка интерполяционного полинома Лагранжа сте- 
пени и-—1 функции (Г) вида [10] 


("—1)! 


запишем минимизационное условие приближения функции /(Х) в виде 


п 1 —1 
С 
ВУ [ОО ат, ри 
А=0 И 


Е 2 
С) +1 Хсыюк-ю- по @ = 


0ЕА=0 


.- 


та "(6 ЧЕ й(пх-К)| > ИМП; 


С: 


1 Ш Ш 
=| р -о = ) 
0 К=0 п 


а А 


[2+ - 2+1 


где й(г) = 2 
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В результате для нахождения С,, К=0,п 


ах(С) 


гебраических уравнений ———^ =0, р=0,п 
р 


Получим, положив С, - ДЖ ЕЕь» п-Д#], вод = [вок К)й(их- руах, 0<^, р<п 


получаем систему линейных ал- 


>в, вв = ор >| еЕ-и о ы 
К=00 К 


п 
Из системы (15) вытекает, что, если Г") =0— =, =0, то С; = Д,, то есть 
для /(х) =Бу +Ьх МНК при указанной Й(Г) является точным. 


№ " 
Если же /(1) =5 +и+5Г, то есть Л (1) =М =соп5 , то Е; „, входящие в 


правые части имеют вид 


0, -р|>2, 
кр =-^| вы. й(пх — К)й(их — руах = р, 
25 п из 51 
м 1 
с рб 
30°" 


Докажем эти формулы. 
При | -Р|>2 получаем зир (их — &) ^зирй(их - р) =@ ‚то есть Е‚›„ =0. 
При к _ Р =1 получаем 


№ ПИ 
Евы =-5. | 6-х - р-1их - р = 
р 


Ре 
М т М ЗИ М 
В [(х р)’ (р+1-пх) =- тия т 
21“ т 21? 5 40 п 
> № 
Аналогично, РЁ’, Е 
40п 


При | Е Р = 0 получаем 


р-+1 


М п 
Реза [по Мик ри рук = 


= 


М 12 М1 
2 п 5 пз 30. 
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Тогда система (15) имеет вид 


п п 
>Врлвь = О Ерь, р=0,п, 
К=0 К=0 


0, р-®>2, 
—_ М 1 
где, в случае |< р<п Е, ры: 
пз 40 
М 1 
——-.—Р-А=0 
‘зо’ 
то есть если р *0,п правые части будут равняться 
ры 1 1 1 
хи» >.Е Кр — е “. “. ее 
а Е 30 п 40 12 
Если о К =0, то 
| 2 1 2 
Во =-У [1-х] Мих-ЮНих- ра =-5 (3-х) мехах = 
| 2 0 п К=0, 2 0 ГД 
р=о 
1 
й 1.21 
Е-^ [5?а-вд?% =- № 2:21 _ _ № 
2. 213 51 60? 
Для случая, когда р=0, А =1, то 
2 
К -- Я (+ мы Кй(их —- р = 
Е 
а 
1 
п 1.1! 
М [1 у вклх — Пика = А [@ — их) Утки 
2? 0 2 2 п? 5 
Тогда если р =0 правая часть будет иметь вид 
о № № № 
Уно = Род + = 5 ея 5. 
0 60" 40п 24п 
Аналогично при р=и получаем 
№ № № 
и Пл Ре == 


4073 6003 243. 


Таким в система алгебраических уравнений может быть записана в виде 


М 
Ву уЕр ==, р=0 
уз рб РЕ 


й=0 
р 
УВ Е == 5р=ЬИ-1> 
К=р-1 
Ш 
У Врл=ь з›Р=П 
К=и-1 
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В условиях вычислительного эксперимента, решая СЛАУ вида (3), получаем 
точные значения матрицы А. Для нахождения эквивалентного возмущения М ис- 


пользуем матрицу А,, полученную с учетом различного количества разрядов ман- 
тиссы. Информация о норме матрицы возмущения позволяет из (13) определить 


=. М 
0з_МЕ 
—#+1 А 
Р Е 
Учитывая, что возмущение вычислительного эксперимента получено только за 
счет учета количества разрядов мантиссы, считаем, что величина (14) является отно- 
сительной ошибкой связанной с округлением. 


и далее относительную ошибку из (14). 


Тогда для функции Х(х) =х? / 2 при п=10, р=10О относительная ошибка 
округления, возникающая при нахождении члена 2 = В а (у) в одномерных опе- 


раторах формулы (2), составляет 00?) при #=10, 00“) при #=15, 00-?) 
при #1 =20. Следует отметить, что для большего п ошибка увеличивается. 


Выводы и перспективы дальнейших исследований 


Выполнение различных научных исследований сегодня невозможно без приме- 
нения ЭВМ для автоматизации сложных вычислительных процессов и преобразова- 
ния большого количества информации. Такие преобразования в конечном счете 
сводятся к выполнению последовательности простейших операций, что и приводит к 
необходимости учитывать ошибку округления на каждом этапе. В данной работе 
предложен подход к оценке ошибки округления, которая возникает в ходе решения 


систем Ви = ЕГ (У) и уВ. = Р. (х) путем вычисления обратной матрицы. 


В дальнейших исследованиях авторы планируют рассмотреть метод оценки 
ошибки округления приближения функций двух переменных через одномерные 
операторы с использованием подходов, описанных в [9]. 
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О.М. Гулуп, О.Т. Лагто5й 


ЕзИтайоп оГКоипатз Еттогу оГ Тио ГапаМез Еипспоп5 


Арргомтайоп оГ Опе-Плтепяюопа[ Орега1ог5 

ТБе агис1е 41$сиз5ез арргоасНез №0 аззеззте гоип4ше етог \Теп арр!уше Ше 
тефо4 оЁ арргохипаНоп о мо уапаМез Юпсйоп$ 4ейпеЯ оп а зубет оЁ тасК$ 
регреп@исч]аг зтале Е Ппез \Веп арргохитайоп орегайог 15 отуеп ш Ше Югт оРа Вооеап 
орегаюг зат гоиэВ Аипепз1опа!| орегаог$ асипе оп Фе опе уапаЫе Рапсйоп [1]. 

Оп Фе аззитрНоп Ша Фе гоип4те еггог оссигз$ ш Фе сасуайоп оР шуегзе 
тайлсез, сопз14еге4 &\м\о угауз оР Нпашсе {ве шуегзе тайлх: ассог4тсе {0 фе Вгзё теоа 
ОЕ зо!ушФ Ппеаг а|сефгалс едиайопз$ зу%етз Ах=ь, ай 15 зеагсв АТ! Базе оп Ше 
десотроз оп ое таблх А шю Шсюгз, Ше зесопа тео 1$ Базе оп Фе аззатриоп 
фа фе тах А! 13 Ме ошу зо шоп ое тайлх едианоп АХ =Е. 

п еасЬ сазе, е шедиаШу олуеп 10 езитае фе гейайуе еггог оЁ зо!уше зузет$ 
Ах =Б Бу шуегзе таблх ассоташе 0 [7], [8]. 

ТБе гезиИ$ оР патенса| ехрегипеп$ Вауе зво\й ФаЁ Фе опе уапа Ме Рапсйоп 
(<) = х2 |2 \Теп и =10, р=10 теайуе гоипаше еггог 15 00-7) аё#=10, 00“) 


аё1=15, 0а0-?) аё 1=20 (Е — Фе патьег о 41215 айег те деспта1 роте (ш Ше сазе 
ог Нхед-роше) ог Ше тапИ$за (ш Фе сазе оР Йоайп-рош) р-агу зубет оЁ апу 
сотрщег сасианоп). ш ад доп, Юг тоге и егтгог шсгеазез. 


Статья поступила в редакцию 20.12.2013. 
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